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Abstract 


We  present  a  variant  of  the  Ab  Intio  Fragment  Orbital  Theory  where  the 
correlation  problem  in  the  excited  states  is  addressed  in  terms  of  the  Coupled  Cluster 
theory  at  the  singles  and  doubles  level .  In  contrast  to  the  earlier  formulations,  the  present 
theory  includes  triple  excitations  in  the  excited  states  and  double  excitations  in  the 
ground  state.  Moreover,  beeause  of  the  inherent  nature  of  Coupled  Cluster  techniques,  the 
eorrelation  thus  included  is  size  consistent.  The  present  theory  is  applied  to  the  study  of 
the  excites  states  of  polyene  chains  of  various  lengths. 


I.  Introduetion 

In  several  previous  reports  [1]  we  have  presented  a  fragmented  ab-initio  approach, 
which  we  have  ealled  the  Ab  Intio  Fragment  Orbital  Theory  (AFOT),  that  consists  of 
treating  a  given  large  molecule  fragment  by  fragment.  The  fragments  are  ehosen  as 
ehemically  viable  units  (such  as  the  functional  groups),  “rectified”  by  adding  hydrogens 
and  then  processed  via  the  Hartree-Fock  methods  to  obtain  the  corresponding  occupied 
and  virtual  orbitals  (called  Fragment  Molecular  Orbitals  —  “FMO’s”).  A  truncated  set  of 
these  orbitals  near  the  top  of  the  valence  band  are  chosen  as  the  active  FMO’s  and  are 
used  to  describe  the  low-lying  excitations.  This  is  accomplished  in  three  steps:  First,  both 
for  representing  the  one-electron  terms  accurately  as  well  as  computing  the  molecular  1- 
and  2-electron  integrals  between  the  nearest  and  also  the  next-nearest  fragments,  new 
Flarteree-Fock  calculations  are  performed  on  what  are  called  segments.  These  are  huilt  by 
surroxmding  each  fragment  by  its  nearest  neighbors.  Secondly,  using  a  Configuration 
Interaction-Singles  (CI-S)  procedure,  the  excited  states  are  generated,  which  are  then 
corrected  by  including  relaxation  terms.  The  method  has  been  found  quite  successful  for 
low-lying  excited  states  (where  the  single  excitations  from  the  Hartree-Fock  gound  state 
account  for  most  of  the  excitation).  They  have  also  been  used  successfully  to  predict 
optical  properties. 

In  extending  the  method  to  organo-metallics  (e.g.  Pt  compoxmds)  [2]  it  was  found, 
however,  that  the  fragment  orbitals  as  obtained  from  the  fragment  calculations  are  not 
adequate  for  acceptable  aeeuracy  for  the  description  of  the  exeited  states.  It  is  necessary 
to  use  instead  their  corresponding  modified  forms  from  their  respective  segment 
calculations.  Clearly  this  new  situation  results  from  a  strong  influence  of  the  neighbors  on 
the  nature  of  the  active  FMO’s.  Since  the  FMO’s  are  now  essentially  derived  from  the 
segment  calculations,  the  corresponding  1-  and  2-  molecular  integrals  can  only  be 
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computed  by  a  second  pass  of  generation  of  the  AO-supermatrices.  This,  however, 
increases  the  AFOT  computation  times  only  slightly,  since  the  bulk  of  the  computing 
time  is  spent  in  performing  the  4-index  transformation  involving  from  AO  to  FMOs’, 
which  step  is  not  altered. 

For  achieving  chemical  accuracy  in  excitation  energies,  we  have  to  consider  two 
types  of  corrections  to  the  CI-S  states  (see  above):  At  the  level  of  the  segment 
calculations,  the  relaxation  corrections  are  to  be  evaluated  for  the  non-excited  occupied 
orbitals.  These  have  been  successfully  carried  out  in  the  past  by  a  second-order 
perturbation  approach,  using  a  truncated  excitation  space.  At  the  global  level,  where  all 
the  active  FMO’s  participate,  the  relaxation  is  obtained  in  terms  of  single  excitations  on 
the  CIS  states.  These  are  in  fact  double  excitations  with  respect  to  the  closed-shell  H.F. 
ground  state.  We  do  not,  however,  allow  these  excitations  in  the  ground  state,  since  these 
excitations  introduce  “correlation”  in  the  Hartree-Fock,  while  they  do  not  on  the  CI-S 
excited  states. 

However,  even  though  in  most  cases,  direct  correlation  terms  do  not  substantially 
alter  the  excitation  energy,  for  low  excitation  energies  they  may  not  be  always  negligible. 
In  order  to  include  the  correlation  terms,  we  have  also  to  ensure  that  they  are  size- 
consistent.  This  is  achieved  most  elegantly  by  treating  the  correlation  via  the  Coupled- 
Cluster  (Singles  and  Doubles)  (CCSD)  theory  [3-5].  In  the  AFOT  framework,  the  most 
important  contribution  of  the  CCSD  treatment  of  the  correlation  takes  place  at  the 
“global”  level  of  calculations  where  all  the  active  FMO’s  are  involved.  We  will  therefore 
restrict  the  application  of  CCSD  only  to  the  correlation  contribution  involving  the  active 
FMOs’. 

We  present  below  a  step-by-step  description  of  the  above  “two-pass”  AFOT 
procedure.  We  also  carry  out  calculations  using  the  present  AFOT  procedure  on  the 
excited  states  of  trans-polyene  chains  of  various  lengths  and  compare  them  with  other 
theoretical  calculations. 


11.  Coupled-Cluster  Wavefunction 
and  Correlation 

The  CCSD  ansatz  for  the  ground  state  is  as  follows: 

T  =  exp(r)|0)  (1) 

where  |0>  is  the  Hartree  Fock  wavefunction  and  T  consists  of  single  and  double 
excitations: 


T  = 


r,+7',  =  Xrx/+ 


ia 


i>  j\a>b 


(2) 
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where  the  ‘t’ s  are  the  excitation  amplitudes. 

These  are  determined  from  the  corresponding  Schroedinger  equation 

//j>P)  =  £ec|>f')  (3) 

with 


^cc 


0Hexp(T)0 


=  E„,+ 
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Starting  from  an  initial  solution 
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these  amplitudes  are  determined  by  the  Galerkin’s  procedure  by  left-multiplying  Eq.  (3) 
by  <0|  i"*"  a  and  <0|(ra)(j^b)  and  solving  the  resulting  equations  iteratively.  The  starting 
solutions,  Eq.  (4),  clearly  represent  the  MP2  approximation  in  an  MBPT  approach. 

In  describing  the  excited  states,  we  shall  follow  the  CCSD  equation  of  motion 
approach  as  proposed  by  Bartlett  and  co workers  with  some  modification.  We  assume  the 
zero-th  order  excited  state  wavefunction  in  the  form: 


M 


(5) 


where  |p>  is  an  excitation  determinant  or  a  spin-symmetrized  combination  of  several 
determinants.  The  corresponding  CCSD  wavefunction  is  given  by 


't'«)  =  exp(r) 


'P. 


(0) 


(6) 


The  total  CCSD  energy  is  given  by: 

E  =  (7) 

where  H  is  given  by: 

H  =  H  +  [H,T]+\lH,T\T]+\^Hj\T\T]  +  i^^H,TlT\T\T]  (8) 

In  the  method  followed  by  Bartlett  et  al  [3,4]  the  t-amplitudes  are  frozen  at  the  values 
obtained  for  the  G.S.  We  would  like  to  depart  from  this  choice  by  arguing  that  while  the 
non-zero  values  of  Ti  result  only  because  of  the  breakdown  of  the  Koopman’s  theorem 
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caused  by  the  presence  of  the  double  excitation  terms  in  the  G.S.,  they  are  largely 
determined  by  the  relaxation  of  the  orbitals  going  from  the  G.S.  to  the  excited  state 
wavefunctions  represented  by  Eq.  (5).  We  solve  for  these  Ti  amplitudes  as  follows: 


ti)  _ 

la 


(9) 


We  shall  take  the  second-order,  T2  coefficients  as 


UJ;ab 


(ia)Ub)\\i^f 


(10) 


The  correlation-relaxation  contribution  will  be  calculated  from  Eq. 


(11) 


Considering  the  fact  that  consists  only  of  single  excitations  from  the  Hartree- 

Fock,  one  can  easily  simplify  the  above  expression  to 


(12) 

AH^,  -2T,HT,\v)  (12a) 


III.  The  Present  Architecture  of  the  AFOT  Procedure 

(a) .The  present  AFOT  code  is  an  add-on  to  the  GAMESS  2004  version  with  the  “control” 
part  of  the  GAMESS  input  retained  intact.  With  the  “AFOT”  keyword  on,  the  program 
reads  in  the  input  for  the  geometry  in  a  slightly  different  format  (see  below). 

(b)  Before  inputting  the  geometry,  the  number  of  active  FMOs’  (Uyai ,  nyir)  are  specified. 
The  geometry  is  given  in  terms  of  fragments  separated  by  a  $  sign.  It  is  implied  that  none 
of  the  fragments  have  common  atoms. 

(c)  The  open-shell  ROHF  wavefunctions  of  the  distinct  atoms  in  their  appropriate  valence 
states  are  obtained  and  the  corresponding  occupied  and  virtuals  orbitals  are  saved  in  the 
RAM. 

(d)  Attach  rectification  hydrogen  atoms  to  the  fragments  to  tie  off  the  dangling  orbitals. 
Form  the  initial  orbitals  by  using  the  information  from  (c).  Obtain  the  RHF  ground  states 
for  each  fragment  and  save  the  resulting  occupied  and  unoccupied  orbitals  in  the  virtual 
memory. 
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(e)  Surround  each  of  the  fragments  by  its  nearest  neighbors,  with  the  latter  rectified  by 
hydrogens.  The  resulting  closed-shell  molecule  is  called  a  segment.  These  segments  are 
now  treated  with  the  RHF  procedure.  The  starting  input  orbitals  are  constructed  by 
including  all  the  occupied  orbitals  of  the  constituent  fragments  except  those  that 
correspond  to  the  orbitals  that  bond  the  rectification  hydrogen  atoms  to  be  eliminated  in 
juxtaposing  the  neighboring  fragments  in  forming  the  segments.  These  are  replaced  by 
what  we  have  termed  the  link  orbitals  obtained  separately  by  diatomic  calculations.  The 
unoccupied  orbitals  of  the  fragments  are  the  stacked  up  on  top  of  the  occupied  ones.  The 
resulting  FMO’s  constitute  the  total  MO-expansion  basis  of  the  segment.  Before  they  are 
used  in  solving  the  segment  Fock  equation,  they  are  symmetrically  orthonormalized  (not 
Schmidt  orthonormalized)  so  they  remain  identified  as  belonging  to  their  original 
fragments  on  equal  footings. 

(f)  Upon  the  convergence  of  the  RFIF  process  and  the  generation  of  the  segment 
canonical  orbitals,  the  active  FMO’s  (AFMO)  localized  in  the  central  fragment  (defined 
as  the  fragment(s)  that  are  not  linked  to  any  rectification  hydrogen  atom)  are  then  derived 
as  follows:  We  expand  the  anticipated  FMO  in  terms  of  the  non-core  orbitals  (both 
occupied  and  virtual)  of  the  central  fragment  in  question  by  a  least  square  fitting  process 
by  minimizing 


L  =  r  (13) 

iciA  idC 

where  {A}  represent  the  space  of  the  active  MO’s  of  the  segment  and  {C}  represent  the 
part  of  the  non-core  orbitals,  associated  with  the  central  fragment  of  interest.  Both  the 
‘a’s  and  the  ‘b’s  are  varied  in  the  least  square-fitting  process.  This  yields  the  optimal 
active  occupied  FMO’s  for  the  central  fragment.  We  obtain  the  excited  FMO’s  as 
follows.  Let  (j)oc  denote  the  highest  occupied  AFMO.  The  corresponding  (j)ex  is  determined 
by  minimizing  the  single-excitation  energy: 

AEex~Sex  “  8oc  “Jex.oc  2Kex,oc  (14) 


This  amounts  to  solving  the  excited  Fock  equation: 

(Fo-.X^2i5C)(l)ex  =  S(j)ex  (15) 

where  (|)ex  is  expanded  in  terms  of  the  unoccupied  space  of  the  central  fragment,  9C 
being  the  Coulomb  and  exchange  operators.  If  it  is  desired  to  have  multiple  excited 
orbitals  <j)ex,  higher  roots  of  the  same  equation  are  selected. 

The  above  AFMO’s  along  with  the  RHF  Fock  matrix  of  the  segment  are  now 

saved. 
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(g)  Upon  completion  of  all  the  segment  calculations,  we  revisit  each  one  of  them  (in  the 
same  sequence)  in  the  “second”  pass  of  the  calculations.  We  recall  the  appropriate  Fock 
matrix  and  all  the  AFMO’s  pertaining  to  all  the  fragments  of  the  segment  in  question. 
Consider  the  total  energy  resulting  from  excitation  of  one  of  the  occupied  AFMO’s  of  the 
central  fragment  to  one  of  the  excited  orbitals  (|)k  belonging  to  any  of  the  fragments 
(central  or  otherwise): 

AEik”  Eq  +  Sk  ~  Si  “Jik"^  2Kik 

The  contribution  from  the  relaxation  of  all  the  occupied  orbitals  except  (t)i  is  given  by  the 
second  order  perturbation  expression: 


AE 


relax  Ak 


-I 


Sk-s, 
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where  j  goes  over  the  occupied  space  excepting  i  and  1  goes  over  different  parts  of  the 
unoccupied  space;  for  j  representing  core  orbitals,  the  whole  imoccupied  space  is  used 
while  for  the  valence  orbitals,  the  unoccupied  space  excludes  all  the  excited  AFMO’s. 


These  relaxation  contributions  are  saved  for  the  future  global  calculations  (see 
below).  We  also  compute  and  save  the  1-  and  2-el  AFMO  integrals,  Fo,pq  and  Ipqrs ,  over 
the  AFMO’s  p,q,r,s. 


(h)  Finally,  the  “global”  calculations,  which  lead  to  the  total  molecular  excited  states  as 
well  as  the  optical  properties,  are  carried  out.  First  we  have  to  evaluate  both  the  2-el 
integrals  as  well  as  the  relaxation  contributions  that  are  not  already  calculated  at  the 
segment  level.  In  this  work,  we  employ  the  simplest  electrostatic  approximation  for  the  2- 
electron  integrals: 


Ipqrs  SpqSrs  /Rpqrs 


where  Rpqrs  is  the  ‘average’  distance  between  the  overlap  pq  and  rs. 


The  relaxation  contributions  for  “distanf  ’  charge  transfers  are  found  by  extrapolation. 

The  first  step  in  global  calculations  is  to  recall  the  saved  1-  and  2-el  integrals  as 
well  as  the  relaxation  contributions  and  carry  out  a  configuration  interaction  (singles) 
(CIS)  calculation.  All  single  excitations  within  the  active  FMO  space  are  included.  Next 
we  compute  the  relaxation  contributions  for  a  CIS-state  of  interest  a  s  given  by 

AErelax  E  Ajk  AErelax,ik  (l^) 
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where  Aik  is  the  Cl-coefficient  for  the  single  excitation  i^k.  This  is  followed  by  a  CCSD 
calculation  yielding  the  contributions  of  the  singles  and  double  excitations  on  the  CIS 
states. 

IV.  Calculations,  Results  and  Discussion 

We  employ  the  above  formalism  to  the  study  of  the  excited  states  (particularly  the 
l^Bu  state)  of  the  trans-polyene  systems  C2nH2n+2,  n=l,10.  The  basis  set  chosen  for  all  of 
our  calculations  reported  in  this  work  is  the  Gaussian  ‘6-31’  set.  In  Table  I  we  summarize 
the  results.  We  present  the  results  for  the  2*Ag  state  (which  lies  consistently  higher  than 
the  Bu  state  because  of  the  absence  of  any  ‘triple-triple’  excitations  in  our  calculations) 
to  indicate  the  narrowing  of  the  2^Ag-l^Bu  gap.  The  correlation  (resulting  from  MP2- 
lecel  double  excitations  on  both  ground  state  and  the  excited  state  (IBu))  seems  to  cancel 
out  routinely  for  all  sizes  of  polyene  systems  considered.  The  contribution  from  the 
‘relaxation  terms’  at  the  global  level  seems  to  be  negligibly  small.  We  must,  however, 
note  that  this  observation  pertains  only  to  the  global  level,  namely  involving  the  homo- 
lumo  orbitals.  There  is  considerable  relaxation  contribution  at  the  segment  level  (  of  the 
order  of  1-2  eV)  from  the  core  orbitals  (see  discussion  above). 
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Table  1.  The  calculated  l^Bu  and  2^Ag  energies  for  various  size  polyenes 


Polyenes 

1  l*Bu  Energy 

1  (eV)  1 

AEmP2  I 
(eV)  1 

AEccsd  I 
(eV)  1 

2  ^Ag  Energy 
!(eV) 

CgHio 

4.46 

<.001 

.08 

6.07 

C12H14 

3.74 

.001 

.11 

4.23 

C16H18 

3.62 

.001 

.12 

4.01 

C20H22 

3.55 

.001 

.12 

3.93 
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